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g ; Abstract 



An upper limit on the Casimir force is found using the dielectric functions 
of perfect crystalline materials which depend only on well defined material 
constants. The force measured with the atomic force microscope is larger 
than this limit at small separations between bodies and the discrepancy is 
significant. The simplest modification of the experiment is proposed allowing 
q-i to make its results more reliable and answer the question if the discrepancy 

has any relation with the existence of a new force. 



The Casimir force between closely spaced macroscopic bodies is the effect 
of quantum electrodynamics (QED) and for that reason could be predicted very 
accurately. In the rigorous Lifshitz theory ||, |3| the force is defined by the optical 
properties of used materials. Knowledge of these properties is the weakest element in 
the theory restricting the accuracy that can be achieved. Though the measurement 
of the Casimir force is not the best way to test QED, such experiments are of great 
importance because they are sensitive to the presence of new fundamental forces Q 
predicted in many modern theories (see, for example, and references therein). 
To distinguish a new force from the background, we should be able to calculate the 
Casimir force with a precision better than the experimental one. In the series of 
recent experiments this force has been measured with the torsion pendulum (TP) 
@ in the range of distances 0.6 — 6 \xm and with the atomic force microscope (AFM) 
|7|, U in the range 0.1 — 0.9 fim. The corresponding precisions were 5% and 1%, 
respectively. 

The force per unit area between parallel plates arising as a result of electromag- 
netic fluctuations at nonzero temperature T is generalized by the Lifshitz theory || , 
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where the plate material is taken into account by its dielectric function at imaginary 
frequencies e 
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Here "prime" means that n = term is taken with the coefficient 1/2, a is the 
distance between bodies and 
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The Casimir result F£ (a) = 7i 2 hc/240a A |JXJ] is reproduced from (|IJ) in the limit 
e — > 00 and T —>■ 0. The function e (i( n ) cannot be measured directly but can be 
expressed via imaginary part of the dielectric function on the real axis with the help 
of the dispersion relation 
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Information on Ime (u) can be extracted from the data on reflectivity and absorp- 
tivity of electromagnetic waves for a given material. 

In the experiments the force is measured between metalized disc and sphere 
because for two plates it is difficult to keep them parallel. For this configuration ([!]) 
has to be modified with the help of the proximity force theorem (PFT) which is 
true for R^> a, where R is the radius of curvature of the spherical surface. Applying 
PFT to (0) one can find the force between sphere and plate as 27tR J F pl (a) da. The 
integration gives 
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This expression differs from those used in || and |T2| in two respects. First, in the 
cited papers the integration connected with the PFT was not done analytically that 
complicated numerical analysis. Second, the zero temperature limit has been taken 
when one can change the sum over n in (|J) by the integral over (. This limit was 



also considered in [13], though the PFT integral was evaluated explicitly. It seems a 
reasonable approximation at small separations because the temperature correction 
is proportional to (kTa/hc) 3 |10[ and is small. However, one should remember 
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that this correction has been found in the limit of ideal conductor e — > oo. For a 
real conductor it can behave as kTa/hc and be important. We have computed the 
force according to and with the integral instead of the sum and found that the 
difference at the smallest distances tested in the AFM experiments exceeds 4 pN in 
contrast with the conservative estimate for the experimental errors 2 pN ||. 

In the AFM experiments an additional Auo^PdoA layer of 20 nm or 8 nm 
H thick was on the top of Al metallization of the bodies to prevent aluminum ox- 
idation. It has to be included into consideration. This layer is transparent for the 
electromagnetic waves with high frequencies ~ c/a since the absorption is propor- 
tional to Ime (u) which is small for u ~ c/a and for this reason the layer was ignored 
in |7|, ||. However, the force depends on e(i() for which the low frequencies dominate 
in (^D because of large Ime (u) and that is why we cannot neglect the Au/Pd layer. 
To take it into account, one has to generalize expression for the force ([I]) to the case 
of layered bodies. Suppose that the top layer has the thickness h and its dielectric 
function is E\. The bottom layer is thick enough to be considered as infinite and let 
its dielectric function be e 2 - The method described in || for deriving Eq.(|l]) can be 
easily generalized for layered plates. We have to add only the matching conditions 
for the Green functions on the layers interface. The result will look exactly as (|]) 
but with more complex Cri )2 : 

Q _ (gi + gg) (p + gi) e^ lfe / c + (gi - s 2 ) (p - gi) e~^ sife / c 
1 ~~ (si + s 2 ) (p - si) e^ h / c + (si - s 2 ) (p + si) e-^ sih / c ' 



_ (gggi + gigg) (giP + gi) e c " Slfe / c + (e 2 si - gigg) (giP ~ gi) e~ c " Slfe / c 
2 ~" (e2Si + 5iS2)(£ip-si)e^ s i /l / c + (e 2 si-5iS2)(giP + si)e^™ s i /l /c' 

where si j2 are defined similar to s in (Q). The force between plate and sphere is 
given by (§]) with the above G\^- Qualitatively the effect of the top layer will be 
negligible if hiO\ p jc <C 1 where uj\ v is the plasma frequency of this layer. For typical 
plasma frequencies ~ 10 16 s~ l it is definitely not the case even for h = 8 nm. The 
force between layered bodies was found also in |Tj| with a little bit different technic 



but it was not used there for actual calculations. 

Now we are able to evaluate the Casimir force in real geometry of the experiments 
if there is information on the dielectric functions of the used materials: Au, Al, and 
AuosPdoA alloy. Strictly speaking, one has to measure these functions in wide range 
of wavelengths on the samples which are used for the force measurement. It was 
not done in all experiments and to draw any conclusion from them we have to make 
some assumptions on e (lu). At low frequencies Au and Al are well described by the 
Drude theory, where 

e(u) = l- . Up . , . (6) 
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Here ui p is the free electron plasma frequency and u T is the Drude damping frequency. 



A simple test for validity of (|g) is the behavior of the material resistivity [14] that 
is defined as 

p{u) = Im— r — 1 = (7) 

where Eo is the free space permittivity. The resistivity is frequency independent 
within the Drude approximation. For crystalline samples of Au and Al (the entries 
2 in Table 1) the frequency behavior of the resistivity and Ime (u) are shown in Fig.l. 



The data on the dielectric functions were taken from [15|], where the data from many 
original works are collected. Palladium definitely cannot be described by (El) at any 
frequency. However, it is known experimentally that amorphous metallic alloys 
such as Au/Pd can be described by the Drude approximation WM. The physical 
explanation for this is associated with large Drude damping of the compounds like 
Au sPd 0A . Eq.(0) allows to use well defined static resistivity p(0) = p instead of 
the damping frequency u T . 

Of course, at higher frequencies when interband transitions are reached the Drude 
approximation fails. Nevertheless, it is very useful since low frequencies dominate 
in the dispersion relation. Extrapolation of @ to high frequencies gives 

The relative error inserted in @ due to extrapolation can be estimated as uj t /ujq, 
where uoq is the frequency of the first resonance for a given metal. The error can be 
as large as 10% but it does not influence significantly in the force. If we will use 
(|) for the force computation and change u p by 5% ( 10% correction to e at all 
frequencies), then the force is changed less than 2%. Moreover, since the interband 
transitions give a correction to (§) which is frequency dependent, it reduces the 
correction to the force further below the experimental uncertainties. The possibility 
to neglect the interband transitions in Al for the force evaluation was noted in 
||. It agrees with our estimate and direct computation using the handbook data 
for Ime{tjj). Therefore, in all cases of interest we can use Eq.(|8]) to describe the 
dielectric function of a material on the imaginary axis. Since the integral in (^) is 
saturated in the low frequency region, we should extract the parameters to p and u T 
from the data for real and imaginary parts of e (u) by fitting them in infrared region 
with (§). 

It is important that optical properties of evaporated (spattered) films can be 
quite different from those of bulk material and depend on the technological details. 
It is known, for example, that the film density is typically 0.7 of that of the bulk 
material if it was not annealed. For the resistivity of spattered and evaporated Au 
films the value po = 8.2 uQ ■ cm has been reported |I| in contrast with the bulk 
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Al 


u p ■ 10" 16 


u T ■ 10" 13 


Po P^ ■ cm 


1* 


1.54 ±0.01 


15.5 ±0.6 


7.39 


2 


2.235 ±0.001 


12.49 ±0.01 


2.83 


3 


2.43 ±0.05 


14.4 ±0.7 


2.76 


4* 


1.63 ±0.03 


18.2 ±0.7 


7.74 


Au 


Up • 10" 16 


u T ■ 10" 13 


Po P^ ■ cm 


1 


1.280 ±0.001 


3.29 ±0.05 


2.27 


2 


1.372 ±0.001 


4.060 ±0.002 


2.44 


3 


1.34 ±0.02 


7.08 ±0.18 


4.46 


4 


1.051 ±0.001 


6.24 ±0.21 


6.40 



Table 1: Drude parameters with statistical errors for Al and Au found by fitting 
e(u) in the infrared range (A > 2 pm). The stars in the first column mark the data 
for film samples. 



resistivity 2.25 pf2 • cm. All these make impossible to use the handbook data for 
reliable calculation of the Casimir force. This conclusion is illustrated by Table [1], 



where the parameters for Al and Au found by fitting the data from [15[ are presented. 

Though we cannot use handbook data to evaluate the force, one can confine it 
for a given experiment. This statement is based on the observation that because 
of better reflectivity the force (f|) increases every time when u p increases or po 
decreases. For us it is important that any technological procedures will reduce u p 
and increase the resistivity po- The perfect crystalline material will have the largest 
plasma frequency and the smallest resistivity and these parameters are well defined. 
The plasma frequency u p is defined by the concentration of free electrons in the 
metal n and their effective mass m* 



e 2 n 



\ m* e e ' 



(9) 



Gold is a good conductor and m* is quite close to the mass of electron. We will find 
the upper limit on the electron concentration if suppose that every Au atom produce 
a free electron. Then for Au plasma frequency one finds Up U = 1.37 ■ 10 16 s -1 . The 
static resistivity can be used to get the damping frequency u T with the help of ([?]) 
at a given u p . For crystalline gold it is p$ u = 2.25 pf2 ■ cm. One can compare 
these parameters with that given in Table 1 to make sure that they correspond 
to the limit values. In the TP experiment || the bodies were covered with Au of 
thickness 0.5 pm that is thick enough to be considered as infinite. Substituting 
the Au parameters into (|j) and calculating the force according to one finds 
the upper limit on the Casimir force in the TP experiment. The residual force 
F exp (aj) — F hm (cii) is shown in Fig.2a, where F exp (aj) are the experimental points 
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at separations eij. The prediction obviously does not contradict to the experiment 
but dealing with the upper limit we cannot conclude that there is an agreement, 
either. 

For the AFM experiments |7J, |8|] the upper limit on the Casimir force is more 
restrictive. The plasma frequency for Al can be restricted using (|9]) if one supposes 
that every atom produces 3 free electrons. It gives Up = 2.40 • 10 16 s -1 that 
coincide with the largest value in Table |l|. The resistivity of perfect crystal is 
Pq 1 = 2.65 pfl ■ cm. Since we successfully predicted the plasma frequencies for the 
best samples of Au and Al, the same way one can use to estimate uj p for Au/Pd. If 
each Au atom gives one and Pd atom gives not more than two free electrons, then 
^Au/Pd _ y gg . io 16 s^ 1 . This alloy is used in microelectronics and resistivity of 

the bulk material is known to be pQ U ^ Pd m 30 /iQ ■ cm [17] in accordance with the 
statement that alloys have large Drude damping. These data allow to find the upper 
limit on the force using with the functions G\ t 2 defined in (|5]). Real surface of the 
bodies is always distorted. The distortion statistics were analysed with atomic force 
microscope B [ill • The force has to be averaged over the distorted surfaces and we 



use for this the procedure developed in |TT[. This procedure seems quite reliable. 



Moreover, the important progress in controlled metal evaporation || allowed to 
reduce the surface roughness to the level when the correction to the force becomes 
practically unimportant for the experiment ||. 

It was indicated J7J that the thickness of Au/Pd layer is less than 20 nm, that 
is why for calculations the conservative value h — 15 nm was chosen. The top layer 
changes the force on 13 pN at the smallest separation. Variation of uip on 10% 
gives only 1 pN change in the force because of screening effect of the top layer. 
The same variation in ujp U l pd changes the force on 2 pN. The resistivity variation 
of the Au/Pd layer on 30% gives 1 pN effect. At larger separations all the effects 
become smaller. All this means that the limit is stable in respect to variation of the 
parameters. It is clear also that the top layer definitely cannot be ignored in the 
force evaluation. The residual force F exp (a) — F hm (a) with the experimental points 
from [J7J is shown in Fig. 2b by triangles. 

In || the assumption of absolute transparency of the Au/Pd layer was not only 
used for theoretical interpretation of the result but also in the procedure of the force 
extraction from the raw data. For this reason we cannot use the points for the force 
directly. Fortunately, it is easy to restore the right data by shifting all the points to 
larger separations on 2h = 16 nm. The result for the residual force F exp (a)—F hm (a) 
with shifted experimental points from || is presented in Fig. 2b by open squares. This 
figure clearly indicates the presence of some unexplained attractive force which is 
decreasing rapidly when the distance between bodies increased. One can speculate 
that the observed discrepancy is explained by a new Yukawa force mediated by 
a light scalar boson but we will not discuss now the restrictions on the Yukawa 
parameters that will be given elsewhere. 



6 



To make the experiment absolutely clear, it is preferable to use Au instead of 
Al metallization because its non-reactive surface has strong advantage over Al. It 
excludes also additional uncertainties connected with Au/Pd layer. One can use 
as well silver or copper but they are not as inert as gold. In practice it is difficult 
to measure the dielectric function at the wavelengths larger than 30 \im but this 
range gives an important contribution to the dispersion relation. That is why the 
material behavior in this range has to be predictable. One can say definitely that 
the materials of platinum group cannot be used since they are not described by (ED 
at low frequencies. An additional advantage of Au metallization is higher density of 
the bodies coating. In this case the hypotetic Yukawa force will be larger roughly 
in (pau/ Pai) 2 ~ 50 times. If the observed discrepancy has relation with the Yukawa 
interaction, the AFM experiment with Au metallization of the bodies will definitely 
reveal this new force even without detailed knowledge of optical properties of the 
metallization. 

In conclusion, we have found the upper limit on the Casimir force that is re- 
alized for perfect crystalline coating of the bodies for which electrical and optical 
properties are well defined. This limit is smaller than the observed force in the 
AFM experiments and the difference far exceeds experimental errors and theoretical 
uncertainties for small separations between bodies. Au metallization of the bodies 
in the AFM experiment will allow to reveal origin of the discrepancy. 
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Figure captions 



Figure 1. Validity of the Drude approximation for Al (triangles) and Au (circles) 
in the infrared range. The resistivity does not depend on frequency (left axis). Solid 
lines (right axis) demonstrate that Ime (u) depends on uj according to @ with the 
parameters given in Table [I] (entries 2). 

Figure 2. The residual force F exp (cij) — F hm (ai) for different experiments: (a) TP 
experiment ||; (b) AFM experiments with the data from (triangles) and from 
P (open squares). 
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